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Q> Abstract 

We review the potential of the LHC to detect a very light CP-odd Higgs boson of the 
NMSSM, a%, through its direct production in association with a bottom-quark pair at 
large tan/3. We also review the LHC discovery potential of the two lightest CP-even 
Higgs states, decaying into two lighter Higgs states or into the lightest CP-odd Higgs 
state and the Z gauge boson. 
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1 Introduction 

The Minimal Supersymmetric Standard Model (MSSM) [I] is probably one of the most 
studied Beyond the SM (BSM) scenarios. However, this model suffers from two critical 
flaws: the /i-problem [2] and the little hierarchy problem. The former flaw results from the 
fact that the Superpotential has a dimensional parameter, \i (the so-called 'Higgs(ino) mass 
parameter'), whose natural value would be either or mp\ (the Planck mass). However, 
phenomenologically, in order to achieve Electro- Weak Symmetry Breaking (EWSB), fi is 
required to take values of the order of the EW scale or possibly up to the TeV range. The 
latter flaw emerged first from LEP, which failed to detect a light CP-even Higgs boson, 
h, thereby imposing severe constraints on m^. For this kind of Higgs state to pass the 
experimental constraints, large higher order corrections from both the SM and SUSY particle 
spectrum are required. The largest contributions come from the third generation, top quarks 
and squarks. However, these required large corrections seem quite unnatural. Recall in fact 
that at tree level the lightest CP-even Higgs boson mass of the MSSM is less than Mz- Even 
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recent LHC results, hinting at the possible existence of a SM-like Higgs state with mass of 
124-126 GeV [21 Sj, weaken the MSSM assumption, as such mass values are really extreme 
in such a SUSY realisation, towards the very end of the allowed mass range. 

The simplest SUSY realisation beyond the MSSM that can solve these two problems at 
once is the NMSSM (for reviews see [HE])- This scenario includes a Higgs singlet Superfield 
in addition to the two MSSM-type Higgs doublets, giving rise to seven Higgs states: three 
CP-even Higgses (mhi < m h 2 < m h 3 ), two CP-odd Higgses 01,2 ijn ai < m a2 ) and a 

pair of charged Higgses h ± . When the scalar component of the singlet Superfield acquires a 
Vacuum Expectation Value (VEV), an 'effective' //-term, /z eff , will be automatically generated 
and can rather naturally have values of order of the EW to TeV scale, as required p 7 ]. In 
addition, in the NMSSM, the little hierarchy problem can be relieved [H [9], since a SM-like 
scalar Higgs boson in the NMSSM context requires less (s)quark corrections than those in 
the MSSM or it can have mass less than the LEP bound due to unconventional decays over 
some regions of the NMSSM parameter spaced). In fact, currently, the NMSSM can also 
explain not only the LHC excess [TUl [EE] a possible LEP excess and is definitely preferred 
by EW global fits [12 H31 E] . 

2 The NMSSM Superpotential 

The Superpotential of the NMSSM is given by 

W = h u QH u U c - h d QH d D c - h e LH d E c + XSH u H d + ^kS 3 , (1) 

where h u , h d , h e , A and k are dimensionless couplings. The term XSH u H d has been intro- 
duced to solve the /z-problem of the MSSM Superpotential. However, the Superpotential in 
Eq. (1) without the term |k5* 3 gives rise to an extra global U(l) symmetry, the so-called 
Peccei-Quinn symmetry U(1)pq [151 US]- Once the Higgs bosons take on VEVs, this sym- 
metry will break spontaneously and lead to the appearance of a CP-odd scalar, called a 
Peccei-Quinn axion. In fact, this axion has not been seen experimentally. In addition, there 
are severe astrophysical and cosmological constraints on A, that is 10 -7 < A < 10~ 10 [17J. 
These constraints necessitate a very large value of < S > in order to solve the /i-problem. 
So, this is not a satisfactory way to solve the latter. 

One elegant way to solve the /z-problem is to break the U(1)pq by introducing an ad- 
ditional term in the Superpotential. This is the last term in Eq. (1) and consequently the 
axion can be avoided. However, introducing this new term in the Superpotential enables 
one to break the PQ symmetry but the Superpotential still have a discrete Z 3 symmetry. 
This discrete symmetry is spontaneously broken when the additional complex scalar field 
acquires a VEV and that will lead to the domain wall problem. That is, during the EW 
phase transition of the early universe, this broken symmetry causes a dramatic change of 
the universe evolution and creates unobserved large anisotropies in the cosmic microwave 
background [IS] . 

3 We will give more explanations of this in Sec. 4. 
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In order to solve the domain wall problem, one needs to break the Z3 symmetry by 
introducing higher order operators at the Plank scale. However, these operators generate 
quadratic tadpoles for the singlet. So, one also needs to impose a new discrete invariance, a 
Z 2 symmetry, on these operators in order to get rid of the dangerous tadpole contributions, 
see [3] for more details. 



3 The Higgs sector of the NMSSM 

The NMSSM Higgs sector contains two Higgs doublets and one Higgs singlet: 

H d = ( 5 ), H u = ( H J Q V S. (2) 



The scalar potential for the Higgs fields can be written as [T9] : 

Vh = V f + V d + Koft, (3) 

where 

V F = \XS\ 2 (\H U \ 2 + \H d \ 2 ) + \XH u H d + kS 2 \ 2 , (4) 

Vn = 9 -^(\H d \ 2 - \H U \ 2 ) 2 + lg 2 \HtH d \ 2 , (5) 

V so{t = m\HlH u + m 2 H H\H d + m\&S + {XA x SH u H d + ^A K S 3 + h.c). (6) 

To generate EWSB, the Higgs fields should have VEVs. In fact, if one assumes that the 
VEVs are real and positive, they can be described by 

At the physical minimum of the scalar potential, Vh, the soft mass parameters of the 
Higgs fields are related to the VEVs through the following relations [T9] : 

2 9 2 i, 



m H u = -g{V d - VJ — -A V d + ~{VZA X + KV S )A 
2 22 1\22,\ \ A V u V d 

m s = -k v s - -A v + KAv u v d + —=\A\ 

The physical Higgs states arise after the Higgs fields acquire VEVs and rotate away the 
Goldstone modes. As a result, the potential can be written as 



V u 1 2 2 
•Vs 7T A V s, 

v d 2 


(8) 
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V H = m 2 h± h + h- + -(A P 2 )M P (p 2 )+\(Si S 2 S 3 )M S \ S, | . (11) 
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The masses of charged Higgs fields, h^, at tree level are 



ml ± =m\ + M^-^\v) 2 , (12) 



where 



Using the minimisation conditions, one can obtain the mass matrices in the scalar and 
pseudoscalar sectors. First, the mass matrix for CP-even Higgs states at tree level has the 
following entries [19] : 

Man = m\+ f Mf - kxv) 2 ] sin 2 2/3, (14) 



M S i2 = ~\( M l - \<M* ) sin4 /3, (15) 
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Msxs = ~\ (m 2 sin2/3 + 2^) (7^) cos2 ^ ( 16 ) 

M S22 = M|cos 2 2/3 + i(At;) 2 sin 2 2/3, (17) 

M m = \ (4^ ff - ^ - ^f^) (18) 

Msss = \mlsm^^f + 4^, +^f^ - h^s^. (19) 

8 ^ s v a 4 

Second, the mass matrix for CP-odd Higgs states at tree level has the following entries [19J: 



Mpu = m% (20) 
M P12 = if m* A Bw2P - Q^) (21) 



2 V A / V^/ieff 

-Mp 22 = ~ f miBin2/3 + 6^) ^sin2/3 - 3 ^f^. (22) 
8 V A / /i eff A 

To a good approximation, at large tan/3 and large m A , the tree level neutral Higgs boson 
masses are given by the following expressions [19J: 

2 3/t/i e fjv4 K 
"V = ^ > (23) 
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l + -(-^)sin 2 2/3), (24) 
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4 LHC phenomenology of the NMSSM Higgs sector 

Because of the existence of a singlet Superfield in the NMSSM, the latter is phenomenologi- 
cally richer than the MSSM. In fact, the NMSSM has seven Higgs states and five neutralinos 
compared to only five Higgs states and four neutralinos in the MSSM. As a consequence, the 
search for Higgs bosons in the context of the NMSSM at present and future colliders is a big 
challenge and more complicated than in the MSSM. 

It was mentioned before that the mass of the lightest CP-even Higgs boson in the MSSM, 
mft, at tree level should be less than Mz- So, large radiative corrections, mainly from top 
and stop loops, are required to pass the LEP lower limit on the Higgs mass. In fact, to 
achieve this we need large stop masses, which only contribute logarithmically in the loop 
corrections. This large discrepancy between top and stop masses causes essentially a fine 
tuning problem [5] (the aforementioned little hierarchy problem). 

As for the NMSSM, the situation is quite different. Assuming CP-conservation in the 
Higgs sector, the upper mass bound for the lightest CP-even Higgs boson at tree level is 
given by 



The last term in this expression can lift up to 10 GeV higher than the corresponding 
mass of the MSSM. So, smaller loop corrections are required to pass the lower bound on the 
SM-like Higgs mass. However, since the higher order corrections are similar to those in the 
MSSM, the upper mass bound reaches 135 - 140 GeV for maximal stop mixing and tan/3 = 2 
[20], [21] . however, this configuration is already excluded in the MSSM by LEP data. Finally, 
notice that the corrections to the lightest CP-even Higgs boson mass are already calculated 
at complete one loop level [221 [221 [21] and also at the dominant two loop level [21]. 

Furthermore, the most interesting property of the NMSSM that can solve the little hierar- 
chy problem of the MSSM comes from the fact that in large areas of the NMSSM parameter 
space Higgs-to-Higgs decays are kinematically open. For instance, the existence of the light- 
est CP-odd Higgs boson a± with mass less than \m hl is quite natural in the NMSSM, see, 
e.g., [25J. In fact, the Branching ratio (Br) for the decay h\ —> e^ai, Br(hi — > aiai), can be 
dominant in large regions of parameter space and as a result the Br(/ti — > bb) is suppressed. 




(27) 
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This unconventional decay channel is so important as it could explain the 2.3a excess ob- 
served at LEP for a Higgs mass, m#, around 100 GeV as shown in figure [TJ The reduced 
coupling in the figure is defined as follows: 

e 2 =(ff?y. (28) 

\9hzzJ 

Here, gf/zz denotes the SM HZZ coupling while guzz the non-standard coupling. As it is 
clear from the plot the excess occurs when the Bt(Hsm —> bb) times £ 2 gives about 20%. 
In the context of the NMSSM, one can explain this excess in two ways. Firstly, a SM- 
like Higgs boson, /i 1)2 , can decay dominantly into a pair of ai's and so the Br(H — > bb) is 
suppressed [T2] [T3] This scenario can relieve the little hierarchy problem but requires 
that m ai < 2mb- (Notice that this mass region is currently highly constrained by ALEPH 
[26] and BaBar [27J data.) In fact, there is also another possibility in the NMSSM that can 
explain the LEP excess due to the fact that the Br (ax — > 77) can be dominant when the a\ 
is highly singlet and again, as a result, the Br(ai — > bb) is suppressed even with m ai > 2rrib 
[2"5] I2"8"] . Secondly, a CP-even Higgs boson, hi, has a reduced coupling with £ < 0.4 [29|, 
due to the mixing between the Higgs singlet and doublets. Notice that neither the SM nor 
the MSSM can explain such modest excess, as they have a ~Bi(H/h — > bb) which is always 
dominant, hence yielding an excess much above the experimental limit. Besides, the NMSSM 
could also explain the recent excess observed at the LHC for a Higgs mass around 125 GeV 

PU EDI ED 122]. 
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Figure 1: Upper limit on the ratio £ from LEP, where the SM Br(H — > bb) and By(H — > 
t + t~) are assumed. Full line represents the observed limit and dashed line represents the 
expected limit. The green band and yellow band are within 68% and 95% probability, 
respectively [6]. 



6 



The discovery of one or more Higgs boson at present or future colliders will open a new 
era in the realm of particle physics. In fact, many efforts have been made to detect such type 
of particles at colliders. In regard to the Higgs sector of the NMSSM, there has been some 
work devoted to explore the detectability of at least one Higgs boson at the LHC and the 
Tevatron. In particular, some efforts have been made to extend the 'No-lose theorem' of the 
MSSM (recall that this states that at least one Higgs boson of the MSSM will be found at 
the LHC via the usual SM-like production and decay channels throughout the entire MSSM 
parameter space [33l EH E5]) to the case of the NMSSM |251I2EIESIEZIEBIE21I3QII3I]. 
By assuming that Higgs-to-Higgs decays are not allowed, it was realised that at least one 
Higgs boson of the NMSSM will be discovered at the LHC. However, this theorem could be 
violated if Higgs-to-SUSY decays are kinematically allowed (e.g., into neutralino or chargino 
pairs, yielding invisible Higgs signals). 

Because of the large number of input parameters of the NMSSM, it is practically very 
difficult to make a continuous scan over all the NMSSM parameter space. The alternative 
way to do the scan is by resorting to benchmark points in parameter space. (For example, 
for benchmark points in the NMSSM parameter space, see Refs. [36] and [12], for the uncon- 
strained and constrained case, respectively.) Either way, one can distinguish between two 
scenarios in which Higgs-to-Higgs decays are either kinematically allowed or not. 

So far, there is no conclusive evidence that the 'No-lose theorem' can be confirmed in the 
context of the NMSSM. In order to establish the theorem for the NMSSM, Higgs-to-Higgs 
decays should be taken into account, in particular the decay h\ — > a\a\. Such a decay can 
in fact be dominant in large regions of the NMSSM parameter space, for instance, for small 
Ak [25], and may not give Higgs signals with sufficient significance at the LHC. However, a 
very light CP-odd Higgs boson, a 1; can be produced in association with chargino pairs [13] 
and in neutralino decays [H] at the LHC. 

The importance of Higgs-to-Higgs decays in the context of the NMSSM has been empha- 
sised over the years in much literature in all the above respects, see, e.g., Refs. [9| Pf5 | H6 | Wf\ . 
Eventually, it was realised that Vector Boson Fusion (VBF)0 could be a viable production 
channel to detect hi 2 — > cnai at the LHC, in which the Higgs pair decays into jjr + r~ J36JHH]. 
Some scope could also be afforded by a 4r signature in both VBF and Higgs-strahlung (off 
gauge bosons) [JS]. The gluon-fusion channel too could be a means of accessing hi — > a\ai 
decays, so long that the two light CP-odd Higgs states decay into four muons [50] or into two 
muons and two taus [STj . Such results were all supported by simulations based on parton 
shower Monte Carlo (MC) programs and some level of detector response. For a recent survey 
of the 'No-lose theorem' in the NMSSM context, see Ref. [52] . 

Besides, there have also been some attempts to distinguish the NMSSM Higgs sector 
from the MSSM one, by affirming a 'More-to-gain theorem' [251 Ell ESI EH EJ E31 El] (that 
is, to recall, to assess whether there exist some areas of the NMSSM parameter space where 
more and/or different Higgs bosons can be discovered at the LHC compared with what is 
expected from the MSSM). Some comparisons between NMSSM and MSSM phenomenology, 
specifically in the Higgs sectors of the two SUSY realisations, can be found in [55J. 

In this paper, we review the LHC discovery potential for the NMSSM Higgs states assum- 

4 Which is dominated by W^+W^-fusion over ZZ-one. 
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ing as production mechanism of these states associated production with bottom-antibottom 
quark pairs. Generally, the heaviest CP-even Higgs, h 3 , the heaviest CP-odd Higgs, a 2 , and 
the h* 1 states have very large masses, above the TeV scale, in particular at large values of 
tan/3, making their discovery at the LHC very difficult. So, we will focus on the LHC dis- 
covery potential through this production mode of the lightest CP-odd Higgs boson, a±, and 
of the lightest two CP-even Higgs states, h\ and hi- 

5 Parameter space scan 

As intimated already, due to the large number of parameters in the NMSSM, it is practically 
not feasible to do a comprehensive scan over all of them. These parameters can however be 
reduced significantly by assuming certain conditions of unification. Here, since the mech- 
anism of SUSY breaking is still unknown, to explore the NMSSM Higgs sector, we have 
performed a general scan in parameter space by fixing the soft SUSY breaking terms at high 
scale to reduce their contributions to the outputs of the parameter scans. Consequently, we 
are left with six independent inputs. Our parameter space is in particular defined through 
the Yukawa couplings A and k, the soft trilinear terms A\ and A K plus tan/3 (the ratio of 
the VEVs of the two Higgs doublets) and /i e fr = A(S) (where (S), recall, is the VEV of the 
Higgs singlet). In our numerical analyses we have taken m^m^) = 4.214 GeV, m^ olc = 1.777 
GeV, mP° lc = 0.1057 GeV and mf' 6 = 171.4 GeV respectively for the running bottom-quark 
mass and the (pole) tau-lepton, muon-lepton and top-quark masses, respectively. 

We have used here the f ortran package NMSSMTools, developed in Refs. [561 |57jl. This 
code computes the masses, couplings and decay widths of all the Higgs bosons of the NMSSM, 
including radiative corrections, in terms of its parameters at the EW scale. NMSSMTools 
also takes into account theoretical as well as experimental constraints from negative Higgs 
searches at LEP [58] and the Tevatron, including the unconventional channels relevant for 
the NMSSM. Notice that the NMSSMTOOLS version used, version 2.3.1, does not include 
the latest LHC constraints [29]. However, as we shall see below, since we keep the SUSY 
mass scales very high and, over the phenomenologically interesting region to this analysis, 
our hi state is not very SM-like, the parameter points tested here are safely beyond current 
LHC limits. 

We have used the code to scan over the six tree level parameters of the NMSSM Higgs 
sector in the following intervals: 



(Notice that our aim is exploring the parameter space which has very low m ai and one way 
to do that is by choosing A K small, in which case its negative values are preferred [19]. Also, 
notice that small A K is preferred to have small fine-tuning |14].) 

Remaining soft terms, contributing at higher order level, which are fixed in the scan 
include: 

5 We have used NMSSMTools_2.3.1. 



A : 0.0001 
/i eff : 100 - 1000 GeV, 



0.7, k : - 0.65, tan/3 : 1.6 - 54, 
A x : -1000 - +1000 GeV, A K :-10 - GeV. 
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• m Q = m t R = ™l R = m~ L = m fR = I TeV, 

• A t = A b = A T = 1.2 TeV, 

• rriq = rriu R = mg R = mj = m~ eR = 1 TeV, 

• M l = M 2 = M 3 = 1.5 TeV. 

As intimated, we have fixed soft term parameters at the TeV scale to minimise their contri- 
butions to parameter space outputs but changing values of some of those parameters such as 
At could decrease or increase the number of successful points emerging from the NMSSM- 
Tools scans but without a significant impact on the m ai distribution. Also, notice that the 
sfermion mass parameters and the SU(2) gaugino mass parameter, M 2 , play crucial roles 
in constraining tan/3. Decreasing values of those parameters allow smaller values of tan/3 to 
pass experimental and theoretical constraints, however, this is a less interesting region of the 
NMSSM parameter space for our analysis, as our Higgs production mode is only relevant 
at large values of tan/3. The effect of heavy gaugino mass parameters on the outputs, in 
particular m ai , would be small except for M 2 through its effect on tan/3. In fact, when 
tan/3 is large, the sfermion masses should be large to avoid the constraints coming from 
the muon anomalous magnetic moment [59]. The dominant Supersymmetric contribution at 
large tan/3 is due to a chargino-sneutrino loop diagram [60]. Also, notice that the chargino 
masses depend strongly on M 2 . 

Guided by the assumptions made in the reference [36], the possible decay channels for 
neutral NMSSM CP-even Higgs boson h, where h = h 12i 3, and neutral CP-odd Higgs boson 
a, where a = ai j2 , are: 

h,a^gg, h,a^fi + fi~, h,a^-T + r~, h,a^-bb, h,a^tt, 
h,a ->■ ss, h,a ->■ cc, h — > W + W~, h ->■ ZZ, 
h,a —> 77, h, a — > Zj, h, a — > Higgses, h,a — >■ sparticles. 

(Notice that the CP-odd Higgses are not allowed to decay into vector boson pairs due to 
CP-conservation.) Also, notice that here 'Higgses' denotes any possible final state involving 
two neutral or two charged Higgs bosons or one Higgs boson and one gauge boson. 

We have performed a random scan over millions of points in the specified parameter 
space. The output of the scan, as mentioned above, contains masses, Br's and couplings of 
the NMSSM Higgses for all the successful points which have passed the various experimental 
and theoretical constraints. 



6 Inclusive event rates 

For the successful data points, we used CalcHEP [61] to calculate the cross sections for 
NMSSM Higgs production^. Some new modules have been implemented for this purpose. 
We focus here on the process 

99 -> bb ai (29) 



6 We adopt herein CTEQ6L [52] as parton distribution functions, with scale Q — VI, the centre-of-mass 
energy at parton level, for all processes computed. 
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i.e., Higgs production in association with a o-quark pair. (The production mode qq — > bb a\ 
is negligible at the LHC with yfs = 14 TeV.) We chose the production mode gg — > bba\ 
because it is the dominant one at large tan/3. The gluon fusion channel is instead burdened 
by huge SM backgrounds and a x does not couple to gauge bosons in Higgs-strahlung and 
Vector Boson Fusion (VBF) processes due to CP-conservation, see [53] . In addition, Higgs 
production in other modes has been studied before, see for example [38]. In fact, Higgs 
production in association with a bb pair has an extra advantage, whereby the associated bb 
pair can be tagged, allowing a useful handle for background rejection. 

In the NMSSM, the a x state is a composition of the usual doublet component of the 
CP-odd MSSM Higgs boson, cxmssm, and the new singlet component, as, coming from the 
singlet Superfield of the NMSSM. This can be written as [H 



ai = gimssm cos 9 A + a s sin 9 A . (30) 

For very small values of Ak, the lightest CP-odd Higgs, a%, is mostly singlet-like with a tiny 
doublet component, i.e., the mixing angle cos 6* a is small, see the top-pane of figure [2] which 
shows the relation between m ai and cos 9 a- The bottom-pane of the figure shows that the 
Br(a! — > 77) can be dominant in some regions of the NMSSM parameter space with the 
possibility of reaching unity when cos 9a ~ 0. 

To a good approximation, m ai can be written in the NMSSM as [1] 

2 KA K /i eff , 2 9 A\fi cS 2 

m ai = -3^— sm 9 A + — T cos 9 A . (31) 



The first term of this expression is dominant at large tan/3. Furthermore, it is clear that a 
combination of all the tree level Higgs sector parameters affects m ai in general. 



7 Photon and tauon signals of very light CP-odd Higgs 
states of the NMSSM at the LHC 

In our attempt to test the two aforementioned theorems, we consider in this section the case 
of the 77 and t + t~ decay channels of a very light CP-odd Higgs boson. The first mode is the 
most important one to detect a CP-even Higgs boson below 130 GeV in the SM and MSSM 
despite the smallness of its branching ratio, of (9(0.001). In addition, this decay mode gives 
a clean signature and can be resolved efficiently at the LHC. The second one is used in the 
MSSM as a search channel of rather heavy CP-even and CP-odd states, in particular at large 
tan/3, and its exploitation has not been proved at very low masses, say, below Mz- 

In the NMSSM, because of the introduction of a complex singlet Superfield, the lightest 
CP-odd Higgs boson, a±, can be a singlet-like state with a tiny doublet component in large 
regions of parameter space. In this section (and also in the next one) we are looking for direct 
production of the a x rather than looking for its traditional production through h 12 decay. We 
examine the discovery potential of the a\ produced in association with a bottom-antibottom 
pair at the LHC through the 77 and r + r~ decay modes. 
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We will show that in the NMSSM there exist regions of its parameter space where one 
can potentially have a dominant di-photon branching ratio of 0(1) for the lightest CP-odd 
Higgs boson with small mass. This possibility emerges in the NMSSM because of the fact 
that such a CP-odd Higgs state has a predominant singlet component and a very weak 
doublet one. As a consequence, all partial decay widths are heavily suppressed as they 
employ only the doublet component, except one: the 77 partial decay width. This comes 
from the fact that the a\x + X~ coupling is not suppressed, as it is generated through the 
XH1H2S Lagrangian term and therefore implies no small mixing. Although the direct decay 
a% — > x + xr is kinematically not allowed, the aforementioned coupling participates in the 
0177 effective coupling [64]. 

Furthermore, we will show that the r-pair decay can be a promising decay mode for 
detecting the a\ state of the NMSSM with very low mass. The detection of such a very 
low mass Higgs state would then unmistakably signal the existence of a non-minimal SUSY 
Higgs sector. 

Figure [3] shows the distribution of the event rates a(gg — > bbai) Br(ai — > 77) and 
cr(gg — > bba\) Br(ai — > t + t~) as functions of m ai and of Br's of the corresponding channel. 
As expected, the inclusive cross section decreases with increasing m ai , see the top panes of 
the figure. It is worth mentioning that the Br(ai — > 77) can be dominant over a sizable 
expanse of the NMSSM parameter space, which originates from tiny widths into all other 
channels due to the dominant singlet nature of a\ as mentioned in Sec. 5.1. However, the 
dominance of Br(ax — > 77) does not correspond to the region that maximises the yield of 
a(gg — > bba\) Br(ai — > 77), as the maximum of the latter occurs for Br's in the region of 
some 10 -5 to 10 -4 , see the bottom-left pane of the figure. Therefore, one can not take full 
advantage of the phenomenon described in the introduction of this section with respect to 
the singlet nature of the a>i state, at the LHC, which couples to 77 through charginos. Thus, 
if ai were highly singlet, it would be difficult for the LHC to discover this particle as the 
doublet component (necessary to enable a large a\bb coupling at production level) would be 
suppressed. The tension between the two components is such that the cross section times 
Br rates are less than 100 fb. 

The outlook for the t + t~ decay mode is much brighter where the corresponding signal 
rates are at nb level for Br(ai — > t + t~) rj 0.1 or even 10 nb for Br(a x — > r + r~) rj 1, 
see the bottom-right pane of figure |3j Also, notice that such large rates naturally hold for 
different values of m 0l , in the allowed interval, but they decrease with increasing m ai (see 
the top-right pane of this figure). 

In the NMSSM, there is a large area of parameter space where one Higgs state can decay 
into two, e.g., hi — > Oiflj: see figure HJ As it is clear from the top-pane of this figure, the 
majority of points generated here have m/ ll > 110 GeV and m ai < 55 GeV, thereby allowing 
the possibility of h\ — > a\a\ decays. Moreover, this decay can be dominant and can reach 
unity as shown in the bottom-pane of the figure. Despite this, such a decay may not give 
Higgs signals with sufficient statistical significance at the LHC (as discussed in previous 
literature). Therefore, we are well motivated to study, in the fortchomin sections, the scope 
of direct production of the ai state in single mode at the LHC, through gg — > bbai, over 



11 



overlapping regions of NMSSM parameter spacd^. 



8 Muon and 6-quark signals of very light CP-odd Higgs 



The di-muon decay mode has an advantage, that it has a clean signature with excellent mass 
resolution. However, the branching ratio is small in most regions of parameter space 

but this decay mode is enhanced for large tan/3. 

Figure shows the correlations between the a\ mass and the di-muon decay rate. One 
can see from this figure that the Br(ai -> /i + /i~) can be of 0(10%), 0(1%) and 0(0.1%) 
or less for the mass intervals 2m M < m ai < 2m T , 2m T < m ai < 2mi> and 2m& < m 0l , 
respectively. The first region of parameter space (m ai < 2m T ) is rather small, the second 
one (2m T < m ai < 2m&) more significant and the third one (2m& < m ai ) is by far the widest 
one. 

Figure [6] illustrates the distribution of the inclusive event rates as a function of the 
Br(ai — > /U + /i~) and of m ai . It is remarkable to notice that the inclusive event rates are 
sizable in all such mass regions. These event rates reach the 10 4 fb level in the two lower mass 
intervals and the 10 3 fb level in the higher mass range and clearly decrease by increasing 
m ai , as expected. Finally, notice that the mass region below the threshold is very 

severely constrained [65] . 




As for 46-quark final states, at large tan/3 values, the cross-section of the a\ produced 
in association with a bottom-antibottom pair followed by the decay a\ — > bb is strongly 
enhanced, in general. However, since the channel is a 4-quark final state, it is plagued by 
very large irreducible and reducible backgrounds. In this section, we examine whether or 
not the production mode gg — > bba\ — > bbbb can be exploited to detect the a± at the LHC. 
In fact, the existence of 6-jets in the final states offers the advantage of 6-tagging, which can 
be exploited to trigger on the signal and enable us to require up to four displaced vertices 
in order to reject light jets. The ensuing 46 signature has already been exploited to detect 
neutral Higgs bosons of the MSSM at the LHC and proved useful, provided that tan/3 is 
large and the collider has good efficiency and purity in tagging 6-quark jets, albeit for the 
case of rather heavy Higgs states (with masses beyond Mz, typically) [661 Ell- 
Figure [7] illustrates the inclusive signal production cross section a(gg — >■ bba\) multiplied 
by the branching fraction Br(di — > bb) as a function of the Br(a x — > bb) and of m ai and 
the plots in figure [S] display instead the correlations between the a x — > bb decay rate and 
the ai mass (top-pane) and between the a x — > bb decay rate and the ai — > 77 decay rate. 
From a close look at the bottom-left pane of figure it is clear that the Br(ai — > bb) is 
dominant for most points in the parameter space, about 90% and above. In addition, by 
looking at the the bottom-right pane of the figure, it is remarkable to notice that also these 
event rates are sizable in most regions of parameter space, topping the 10 7 fb level for small 
values of m ai and decreasing rapidly with increasing m ai . One can also notice that there 

7 A partonic signal-to-background (S/B) analysis for 77 and t + t~ final states has been done in [25] , 
where extraction of the latter signature was proven for several benchmark scenarios. 



states of the NMSSM at the LHC 
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are some points in the parameter space with m ai between 40 and 120 GeV, as shown in 
the top-pane of figure El in which the Br(ai — > bb) is suppressed due to the enhancement 
of the Br(ai — > 77) (see the bottom-pane of the same figure), a phenomenon peculiar to 
the NMSSM that depends upon the amount of Higgs singlet-doublet mixing^], see [22] • (We 
discussed this in the previous section^.) 

9 The 'No-lose theorem' for NMSSM Higgs discovery 
at the LHC in difficult scenarios 

In this section, we continue to investigate whether or not the 'No-lose theorem' of the NMSSM 
at the LHC can be proven considering Higgs boson production in association with a o-quark 
pair. We do so based on the interesting results obtained in the previous two sections. In this 
section we will, however, no longer consider direct a>i production, i.e., gg — > bba\. Rather, we 
will initially produce either a hi or h 2 , eventually decaying to one or more ai's. In this case 
one may wonder though whether also Higgs boson production in association with a t-quark 
pair could play a role, owning to different couplings of the h% t 2 state to fermions, with respect 
to the cti field. Production rates for gg — > tth\^ were studied in [38], where they were found 
to be very subleading over the entire NMSSM parameter space. 

We will be looking at inclusive event rates in presence of various Higgs-to-Higgs decays, 
namely, — > aiOi, h 2 — > h\h\ and — > Za\. We will also be studying the decay patterns 
of the lightest Higgs boson pairs, a\d\ or h\h\, and of the gauge boson and a light CP-odd 
Higgs boson, Za±, into different final states. Further details on material contained in this 
section can be found in [40 J and [4"T] . 

We have again used NMSSMTools to perform a random scan over the usual parameter 
space, mentioned in section 5, and further required that < 300 GeV, as corresponding 
production rates become negligible for heavier masses. We used CalcHEP [61] to determine 
the cross sections for NMSSM hi t2 production for the following two processes: 

gg — > bb hi and gg —> bb h 2 , (32) 

which were computed separately (i.e., without the interferences emerging whenever hi and 
h 2 have the same decay products). 

The lightest two CP-even neutral Higgs boson masses are given by Eq.(25). Recall though 
that the equation is at tree level, mainly for guidance in interpreting the upcoming figures, 
while NMSSMTools includes radiative corrections as well. 

Figure [9] shows the correlations between all three Higgs masses m ai , m/ ll and m^ 2 . Since 
the successful points emerging from the scan have small values of A, n and also A K , only 
rather small values of m ai are allowed. It is remarkable that the smaller m ai , the smaller 

8 Notice that constraints coming from Tevatron [55J do not affect our results since the singlet field plays 
a primary role in the NMSSM. But under severe conditions such as A — > and k — > 0, the NMSSM and 
MSSM become similar and those constraints may be applied. 

9 A partonic signal-to-background (S/B) analysis for and bb final states has been done in [251 155] , 

extracting both signals for several benchmark scenarios. 
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and rrih 2 too (two top-panes). In the bottom-pane of the same figure, for rrih 2 around 
120 GeV, m/ ll can have values from just above up to slightly less than 120 GeV, showing 
the possibility that the two Higgs states can have the same mass, i.e., m hl ~ m h2 . Notice 
also that the majority of points have between 115 and 120 GeV. 

9.1 Production of h\ and hi decaying into two lighter Higgs bosons 

The production times decay rates of h x and h%, in which hi decays into two lighter ai's and 
hi decays into either a pair of ai's or a pair of /ii's, are shown in figure [TDJ This figure 
displays all the correlations between the three discussed production and decay processes. It 
is quite remarkable that the overall trend, despite an obvious spread also in the horizontal 
and vertical directions, is such that when one channel grows in event yield there is also 
another one which also does, hence opening up the possibility of the simultaneous discovery 
of several Higgs states of the NMSSM (three neutral Higgses at the same time: hi, hi and 
di), an exciting prospect in order to distinguish the NMSSM Higgs sector from the MSSM 
one (in fact, a clear manifestation of a possible More-to-gain theorem being established j^i 

9.2 Production of h\ and hi decaying into a gauge boson and a 
light CP-odd Higgs 

In this subsection, we examine the LHC discovery potential of the lightest two CP-even Higgs 
states hifli followed by the decay h x ^ ~ > Za x . Figure [TT] shows that the production rate for 
the hi, produced in association with a bb pair, is small, topping the 0.001 fb level. Such 
a production rate is presumably not enough to discover the hi at the LHC. The top-pane 
of the figure shows that there is a linear relation between the hi production rate and the 
Br(/ix Za x ) because the production rate u{gg — > bbhi) is nearly constant in our parameter 
space, which has large tan/3. The bottom-pane of the figure shows that the points passing 
the constraints have > 100 GeV. 

Figure [12] illustrates the inclusive hi production rates ending up with Zai — > fi + [i~bb, 
Za x — > fi + fi~T + r~ and Zai H T+r ~ (where j = jet). It is clear that the production and 
decay rates are definitely too small, topping 10 -4 fb for the first and last channels and 10~ 5 
fb for the second one. Such rates are obviously not enough to discover the hi neither at the 
LHC nor at the SLHC with 1000 fb' 1 of luminosity. 

In contrast, the situation for hi is promising as one can notice that u{gg — >• bbh 2 )Br(h 2 — > 
Zai) is sizable, topping the 10000 fb level (figure ITS]) . The highest values of the cross section 
are accompanied by an intriguingly large Br(hi — > Zai), reaching up to 10%. It is clear from 
the top-pane of the figure that the distribution over the branching ratio for h 2 is not as 
uniform as that for the hi because the production rate a(gg —> bbh 2 ) depends strongly on 
the tree level parameters unlike that for h x . The bottom-pane of the figure shows that the 
highest cross section occurs for mh 2 > 220 GeV. 

10 A partonic signal-to-background (S/B) analysis for several a\a\ and h\h\ decays is currently being done 
in [BS]. 
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In order to study the detectability of h 2 decaying into a gauge boson and a light CP-odd 
Higgs state at the LHC, we have calculated the inclusive production rates ending up with 
fi + fi~bb, [i + {i~t + t~ and jjr + r~ (figure [Tj]) . The event rates for these processes are at the 
(9(100) fb level at the most. While clearly this number is not very large, signal events 
may still be detectable at planned LHC luminosities, especially if the background can be 
successfully reduced to manag eable levelH In short, there is a small but well defined region 
of the NMSSM parameter space where the h 2 and a\ states, both with a mixed singlet and 
doublet nature, could potentially be detected at the LHC if 220 GeV < mya < 300 GeV 
and 15 GeV < m al < 60 GeV, in the h 2 — > Zai — > n + ^~bb, h 2 — > Zai — > fi + fi~r + T~ and 
h 2 — > Zai — > jjr + r~ modes, when the CP-even Higgs state is produced in association with 
a bb pair for rather large tan/3. 

10 Conclusions 

The NMSSM has a singlet Superfield in addition to the usual Higgs doublets of the MSSM. 
This singlet gives rise to a more varied phenomenology in the case of the NMSSM, compared 
to that of the MSSM. For instance, this singlet Superfield mixes with the neutral components 
of the doublets, giving rise to one CP-even Higgs, one CP-odd Higgs and one extra neutralino 
in addition to the usual spectrum of the MSSM. Therefore, in the NMSSM, by assuming 
CP-conservation, there are seven Higgses: three CP-even, two CP-odd and a pair of charged 
Higgses. We have investigated whether or not at least one Higgs boson of the NMSSM can 
be discovered at the LHC ('No-lose theorem') and/or is possible to find some regions in the 
parameter space where more and/or different Higgs states of the NMSSM are detectable at 
the LHC, compared to those available within the MSSM ('More-to-gain theorem'). 

Because of the mixing between the Higgs singlet and doublets, Higgs-to-Higgs decays 
are kinematically possible for large regions of the NMSSM parameter space even for small 
masses of the Higgs states, which is impossible in the MSSM. For instance, a SM-like Higgs 
can decay into a pair of the lightest NMSSM CP-odd Higgses. This decay can be dominant 
in sizable areas of the NMSSM parameter space. Such a decay has a significant meaning if 
one notices that it can explain a 2.3a excess occurred at LEP for the process e + e~ — > Zbb for 
M fe ^~ 98 GeV and the 2.6a excess recently emerged at the LHC (primarily in the 77 decay 
mode). Moreover, a SM-like Higgs with mass of order 100 GeV, which has no-fine tunning, 
can naturally occur in the NMSSM and this scenario is preferred by precision EW data. In 
addition, the NMSSM can solve both the /i-problem and the little hierarchy problem of the 
MSSM. 

In the context of the NMSSM, we have proven that a very light CP-odd Higgs state 
with mass m ai < Mz, which has a large singlet component and a small doublet one, can 
be discovered at the LHC via Higgs production in association with a bottom-antibottom 
pair. This mode is dominant at large tan/3. After performing several analyses for signals 
and dominant backgrounds, not documented here yet referred to, we have proven that this 

11 A partonic signal-to-background (S/B) analysis for jjr + r~ final state has been done in [41], showing 
very promising results. 
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production mode is the ideal one to discover the a\ through the following signatures: (i) t + t~ 
decay mode, in which a\ can be discovered with mass up to Mz\ {%%) /U + [i~ decay mode, if 
10 < m ai < 60 GeV. Further, despite the fact that the bb decay mode is dominant in most 
regions of parameter space that have light ai, this channel has huge QCD background and a 
smaller signal-to-background ratio. Finally, we also looked at the detectability of a\ through 
the 77 decay mode but this proved unuseful despite the fact that this decay mode can be 
dominant in some areas of the NMSSM parameter space. 

We believe that the results presented in sections 7 and 8 have a twofold relevance. Firstly, 
they support the 'No-lose theorem' by looking for direct ai production rather than looking 
for its production through the decays h 12 — > OiOi, which may not give a sufficient signal 
significance. Secondly, they corroborate the 'More-to-gain theorem' as such very light a±s 
(with m ai < Mz) are not at all possible in the MSSM. Altogether, the existence of such a 
light neutral Higgs state is a direct evidence for the non-minimal nature of the SUSY Higgs 
sector. 

Finally, we have mentioned in section 9 the importance of Higgs-to-Higgs decays in the 
NMSSM, here occurring after Higgs boson production in association with bb pairs (unlike in 
most previous literature), and have shown that such decays should be taken seriously before 
proving, or otherwise, the 'No-lose theorem'. In fact, we also have shown that such decays 
are dominant in sizable regions of the NMSSM parameter space. We have studied the LHC 
discovery potential of a CP-even Higgs boson hi or h 2 , decaying into a pair of light CP-odd 
Higgses ax's, and also h 2 decaying into a pair of his. We have found that these channels can 
give sizable signal rates, which could allow one to detect simultaneously two Higgs bosons: 
hi and ax, h 2 and ai or h 2 and h\. In addition, we have shown that the LHC has the potential 
to discover the three neutral Higgs bosons at the same time. Furthermore, we have studied 
the LHC discovery potential for hi and h 2 decaying into Zai and have shown that, while 
the discovery of the hi through this channel is impossible, there is a small but well defined 
region of the NMSSM parameter space where the h 2 state could potentially be discovered. 
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Figure 2: The lightest CP-odd Higgs mass m ai and the Br(ai — > 77) plotted against the 
mixing angle in the CP-odd Higgs sector cos#a- 
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Figure 3: The rates for cr(gg — >■ bbai) Br(ai — > 77) (left) and for a(gg —> bbai) Br(ai 
+ r _ ) (right) as functions of m ai and of the Br of the corresponding channel. 
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Figure 4: The lightest CP-odd Higgs mass m ai plotted against the lightest CP-even Higgs 
mass m hl and against Br(/ii — > aiai). 
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Figure 5: The CP-odd Higgs mass m 0l as a function of the Br(ai — > ). 
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$ure 8: The Br(ai — > bb) as a function of the CP-odd Higgs mass m ai and of the Br(ai — > 



25 



140 



120 - 



100 - 




140 



120 - 



100 - 



0) 80 - 

a 




100 120 140 160 180 200 220 240 260 280 

m h2 [GeV] 




20 - 



1 1 1 1 1 1 1 1 1 1 1 

100 120 140 160 180 200 220 240 260 280 300 

m h2 [GeV] 



Figure 9: The correlations between the lightest CP-odd Higgs mass m ai and the lightest two 
CP-even Higgs masses m hl and m h2 and between the latter two. 
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Figure 10: The rates for a(gg — >■ Br(/i 1 — >■ aiOi) versus cr(gg — >■ &&/i 2 ) Br(/t 2 — > 
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Figure 11: The signal rate for a(gg — > bbh\) Br(/ii — > Za\) as a function of the Br(/ii — > Za\] 
and of mh x - 
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Figure 12: The signal rate for a(gg — > bbhi) Bi(hi — > Za\) times ~Bi{Za\ — > bb), times 
Br(Zai — > n + h~t + t~) and times Br(Zai — > jjr + r~) as functions of m hl . 
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Figure 13: The signal rate for a(gg — > 66/12) Br(/i2 — > Za\) as a function of the Br(/i2 — > Za\ t 
and of rrih 2 . 
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Figure 14: The signal rate for a(gg — > bbh2) Br(/i 2 — > Za\) times ~Bi{Za\ — > bb), times 
Br(Zai — > n + h~t + t~) and times Br(Zai — > jjr + r~) as functions of m h2 . 
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